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The ground state of bosonic atoms in a trap has been shown experimentally to display Bose-Einstein 
condensation (BEC). We prove this fact theoretically for bosons with two-body repulsive interaction 
potentials in the dilute limit, starting from the basic Schrodinger equation; the condensation is 100% 
into the state that minimizes the Gross-Pitaevskii energy functional. This is the first rigorous proof 
of BEC in a physically realistic, continuum model. 
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It is gratifying to see the experimental realization, in 
traps, of the long-predicted Bose-Einstein condensation 
(BEC) of gases. From the theoretical point of view, 
however, a rigorous demonstration of this phenomenon 
- starting from the many-body Hamiltonian of interact¬ 
ing particles - has not yet been achieved. In this letter we 
provide such a rigorous justification for the ground state 
of 2D or 3D bosons in a trap with repulsive pair poten¬ 
tials, and in the well-defined limit (described below) in 
which the Gross-Pitaevskii (GP) formula is applicable. 
It is the first proof of BEC for interacting particles in 
a continuum (as distinct from lattice) model and in a 
physically realistic situation. 

The difficulty of the problem comes from the fact that 
BEC is not a consequence of energy considerations alone. 
The correctness JlJ of Bogolubov’s formula for the ground 
state energy per particle, eo (p), of bosons at low density 
p , namely eo (p) = 2ttTi 2 pa /to (with m = particle mass 
and a = scattering length of the pair potential) shows 
only that ‘condensation’ exists on local length scales. 
The same is true || in 2D, with Schick’s formula || 
e 0 (p) = 2irTi 2 p/{m\ ln(pa 2 )|). Although it is convenient 
to assume BEC in the derivation of eo (/?), these formulas 
for eo (p) do not prove BEC. Indeed, in ID the assump¬ 
tion of BEC leads to a correct formula ||] for eo (p), but 
there is, presumably, no BEC in ID ground states J|] . 

The results just mentioned are for homogeneous gases 
in the thermodynamic limit. For traps, the GP formula 
is exact HfiJ in the limit, and one expects BEC into the 
GP function (instead of into the constant, or zero mo¬ 
mentum, function appropriate for the homogeneous gas). 
This is proved in Theorem 1. In the homogeneous case 
the BEC is not 100%, even in the ground state. There 
is always some ‘depletion’. In contrast, BEC in the GP 
limit is 100% because the IV —> oo limit is different. 

In the homogeneous case one fixes a > 0 and takes 
N —► cio with p = IV/volume fixed. For the GP limit one 
fixes the external trap potential U(r) and fixes Na , the 
‘effective coupling constant’, as N —> oo. A particular, 
academic example of the trap is V(r) = 0 for r inside 
a unit cube and V(r) = oo otherwise. By scaling, one 
can relate this special case to the homogeneous case and 


thereby compare the two limits; one sees that the homo¬ 
geneous case corresponds, mathematically, to the trap 
case with this special V, but with Na 3 = pa 3 fixed as 
IV —> oo. Thus, BEC in the trap case is the easier of the 
two, reflecting the incompleteness of BEC in the homo¬ 
geneous case. The lack of depletion in the GP limit is 
consistent with pa 3 —> 0 and with Bogolubov theory. 

We now describe the setting more precisely. We con¬ 
centrate on the 3D case, and comment on the generaliza¬ 
tion to 2D at the end of this letter. The Hamiltonian for 
N identical bosons in a trap potential V, interacting via 
a pair potential v, is 

N 

H = ^(- A i-t-T( r i))-(- X v ( r i~ r j)- W 

7=1 1 

It acts on symmetric functions of N variables ry E R 3 . 
Units in which fi 2 /2m = 1 are used. We assume the 
trap potential V to be a locally bounded function, that 
tends to infinity as |r| —» oo. The interaction potential 
v is assumed to be nonnegative, spherically symmetric, 
and have a finite scattering length a. (For the definition 
of scattering length, see I- 1 or Q.) Note that we 
do not demand v to be locally integrable; it is allowed 
to have a hard core, which forces the wave functions to 
vanish whenever two particles are close together. In the 
following, we want to let a vary with IV, and we do this 
by scaling, i.e., we write v(r) = V\(r/a)/a 2 , where v\ has 
scattering length 1 , and keep v\ fixed when varying a. 

The Gross-Pitaevskii functional is given by 

£ GP [4>] = J (|V</(r )| 2 + U(r)|<^(r )| 2 + < 7 | 0 (r) | 4 ) d 3 r. 

The parameter g is related to the scattering length of the 
interaction potential appearing in (jl|) via 

g = 47rIVa. (2) 

We denote by <fi GF the minimizer of £ GP under the nor¬ 
malization condition f \(f>\ 2 = 1. Existence, uniqueness, 
and some regularity properties of </ GP were proved in the 
appendix of |(J. In particular, (f> GF is continuously differ¬ 
entiable and strictly positive. Of course </> GP depends on 


1 


g , but we omit this dependence for simplicity of notation. 
For use later, we define the projector 

P G P = |^GP )W) GP| ( 3 ) 

It was shown in |Gj] (see also Theorem 2 below) that, for 
each fixed g , the minimization of the GP functional cor¬ 
rectly reproduces the large N asymptotics of the ground 
state energy and density of H - but no assertion about 
BEC in this limit was made in |]]. 

BEC in H/, the (nonnegative and normalized) ground 
state of H, refers to the reduced one-particle density ma¬ 
trix 

7 (r,r ') = N J tf(r,X)tt(r',X)dX , 


where X = (r 2 ,..., rjy) and dX = ]~f ^ 2 d 3 r r 

Complete (or 100%) BEC is defined to be the property 
that ^7 becomes a simple product /(r)/(r') as N —> oo, 
in which case / is called the condensate wave function. 
In the GP limit, i.e., TV —► oo with g = AnNa fixed, we 
can show that this is the case, and the condensate wave 
function is, in fact, the GP minimizer <p GP . 

THEOREM 1 (Bose-Einstein Condensation). For 
each fixed g 

lim -^ 7 (r,r') = 0 GP (r) 0 GP (r') . 

N —>oo iv 

Convergence is in the senses that Trace | -^7 — P GP | —> 0 

and f (^ 7 (r, r') - </> GP (r)</> GP (r ')) 2 d 3 rd 3 r' —> 0 . 

We remark that Theorem 1 implies that there is 100% 
condensation for all n-particle reduced density matrices 
of T, i.e., they converge to the one-dimensional projector 
onto the corresponding n-fold product of cf GP . To see 
this, let a *, a denote the boson creation and annihilation 
operators for the state 4> GP , and observe that 


1 > N~ n (^\(a*) n a n \^) « iV- n (4'|(a*a) n |'F} 

> TV“ n (4'|a*a|1 , ) n -> 1 , 

where the terms coming from the commutators [a, a*] = 1 
can be neglected since they are of lower order as N —> 00 . 
The last inequality follows from convexity. 

Another corollary, important for the interpretation of 
experiments, concerns the momentum distribution of the 
ground state. 

COROLLARY 1 (Convergence of Momentum Distri¬ 
bution). Let p(k) = f y(r, r') exp[?'k • (r — r , )]d 3 rd 3 r / 
denote the one-particle momentum density of At. Then, 
for each fixed g , 

lim T7P( k ) = l? GP ( k )| 2 

N —>oo J\J 


the sense that f -^-p( k) — |^> GP (k) 


d 3 k —* 0. Here, 


(f> GP denotes the Fourier transform of <f> GP . 


Proof. If T denotes the (unitary) operator ‘Fourier 
transform’ and if ip is an arbitrary bounded function with 
bound H^Hoo, then 




Trace [P f ( 7 /N - P GP )^] 


< 11 (f 1100 Trace ( 7 /^ - P GP |, 


whence f p/N — \(f> 


GP12 


< Trace (7 /IV — P GP |. 


QED 


Before proving Theorem 1, let us state some prior re¬ 
sults on which we shall build. Then we shall outline the 
proof and formulate two lemmas, which will allow us to 
prove Theorem 1. We conclude with the proof itself. 

Denote by E CiM (N, a) the ground state energy of H 
and by E GP (g) the lowest energy of £ GP with f |</>| 2 = 1. 
The following Theorem 2 can be deduced from |j;]. 

THEOREM 2 (Asymptotics of Energy Components). 
Let p{ r) = 7 (r, r) denote the density of the ground state 
of H. For fixed g = 47 r./Va, 


lim l E ^ M (N,a)=E GP (g) (4a) 

N—>oo A 


and 

lim lp(r) = |^ GP (r )| 2 (4b) 

iv—»00 iv 


in the same sense as in Corollary 1. Moreover, if 
denotes the solution to the scattering equation for v\ 
(under the boundary condition limiri^oo <pi(r) = 1) and 
s = f |V</ 5 1 | 2 / 47 r, then 0 < s < 1 and 


lim / \S/ ri ^(r l ,X)\ 2 d 3 r 1 dX 

N-> 00 J 

= J jV</> GP (r)| 2 d 3 r + gs J |</i) GP (r)| 4 d 3 r, (5a) 

lim /p(n)|^| 2 d 3 ndX= /p(r)|^. GP (r)| 2 d 3 r, (5b) 

N-> 00 J J 

N 

lim i y 

iV—>00 1 ^ 

3 =2 

= (1 - s)g J |</> GP (r)| 4 d 3 r. (5c) 

Only (Q) was proved in Jg| , but (® follows, as noted in 
Q, by multiplying V and v by parameters and computing 
the variation of the energy with respect to them. 

(Technical note: The convergence in (^) was shown 
in !| to be in the weak L 4 (R 3 ) sense, but our result 
here implies strong convergence, in fact. The proof in 
Corollary 1, together with Theorem 1 itself, implies this.) 

Outline of Proof: There are two essential ingredients 
in our proof of Theorem 1. The first is a proof that 
the part of the kinetic energy that is associated with the 
interaction v (namely, the second term in (^a|)) is mostly 
located in small balls surrounding each particle. More 


J v (ri - rj)|’F(ri,X)| 2 d' : 


ri dX 
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precisely, these balls can be taken to have radius N~ 'Z 17 , 
which is much smaller than the mean-particle spacing 
N- 1 / 3 . This allows us to conclude that the function of r 
defined for each fixed value of X by 

/xW = ^yf <r ’ x)ao (6) 


Writing / x (r) = n fc > 2 0 GP (r fc ).F(r, X) and using that 
F is symmetric in the particle coordinates, we see that 
(ip.j) is equivalent to 

-~ 9 J l^ GP l 4 - °( 1 ) ; ( 10 ) 

where Q is the quadratic form 


has the property that V r /x(r) is almost zero outside the 
small balls centered at points of X. 

The complement of the small balls has a large volume 
but it can be a weird set; it need not even be connected. 
Therefore, the smallness of V r /x(r) in this set does not 
guarantee that /x( r) is nearly constant (in r), or even 
that it is continuous. We need /x(r) to be nearly con¬ 
stant in order to conclude BEC. What saves the day is 
the knowledge that the total kinetic energy of /x( r ) (in¬ 
cluding the balls) is not huge. The result that allows us 
to combine these two pieces of information in order to 
deduce the almost constancy of /x(r) is the generalized 
Poincare inequality in Lemma 2. (End of Outline.) 

Using the results of Theorem 2, partial integration and 
the GP equation (i.e., the variational equation for <p GF , 
see || , Eq. (2.4)) we see that 


J |</> GP (r)| 2 |V r / x | 2 d 3 r dX = gs J \f> GP \ 4 d 3 v. (7) 


The following Lemma shows that to leading order all the 
energy in (Q) is concentrated in small balls. 

LEMMA 1 (Localization of Energy). For fixed X let 


fix = 


jr G R 


3 


min |r — 
fc> 2 



( 8 ) 


Then 


lim 

iV—»oo 


dX f d 3 r|<p GP (r)| 2 |V r /x(r)| 2 =0. 
J fix 


Proof. We shall show that, as N —> oo, 



- 2 g JdxJ d 3 r|<(» GP (r)| 4 |/x(r)| 2 
> ~9 J |<(> GP (r)| 4 d 3 r - o(l) , (9) 

which implies the assertion of the Lemma by virtue of ( 0 ) 
and the results of Theorem 2. Here, fl^ is the comple¬ 
ment of fix- The proof of ( 0 ) is actually just a detailed 
examination of the lower bounds to the energy derived 
in H and [ij, and we use the methods in jjj.jl] . just de¬ 
scribing the differences from the case considered here. 


N 


N 


Q(F) = E / l V ^| 2 n l</> GP (r fe )| 2 rf 3 r fe 

i=l -' n ? fe=1 

N 

+ E / _ r i)l F | 2 n l^ GP (r fc )| 2 <i 3 r fc 

1 <i<j<N J k =1 

N r N 

-2gJ2 / l^ GP ( r i)\ 2 \F\ 2 LI \t GP (r k )\ 2 d 3 r k , (11) 

A 1 J _ 1 


k =1 


with fi? = {(ri,X) G R 3N \ min*^ |rj - r fc | < 7 / 17 }. 

While (]l^) is not true for all conceivable F’s satisfying 
the condition J |E| 2 ]~[^ =1 |0 GP (r^)| 2 cf 3 rfc = 1, it is true 
for an F, such as ours, that has bounded kinetic energy 
(0). Eqs. 14.1 1)—(4. 12), (4.23)-(4.25), proved in Jg| , are 
similar to ([!(]), ( |ll| ) and almost establish (|l^), but there 
are two differences which we now explain. 

(i) In our case, the kinetic energy of particle i is re¬ 
stricted to the subset of R 3W in which min^^j |r» — r^| < 
N~ 7 /i7. However, from the proof of the lower bound to 
the ground state energy of a homogeneous Bose gas de¬ 
rived in [jlj (especially Lemma 1 and Eq. (26) there), 
which enters the calculations in we see that only 
this part of the kinetic energy enters the proof of the 
lower bound — except for some additional piece with 
a relative magnitude e = 0(iV -2 / 17 ). In the notation 
of jl]] the radius of the balls used in the application 
of Lemma 1 is chosen to be R = aE -5 / 17 , which, in 
the GP regime, is R = 0(N~ 7 / 17 ) since, for fixed Na , 
Y = 0(a 3 N) = 0(N~ 2 ). (See Q for a fuller discus¬ 
sion about the choice of R.) The a-priori knowledge that 
the total kinetic energy is bounded by (jTI) tells us that 
the ‘additional piece’, which is e times the total kinetic 
energy, is truly O(e) and goes to zero as N —> oo. 

(ii) In H all integrals were restricted to some arbi¬ 
trarily big, but finite box of size R'. However, the dif¬ 
ference in the energy is easily estimated to be smaller 
than 2 gN x max| r |> ff |0 GP (r) | 2 , which, divided by N, 
is arbitrarily small, since (f> GF ( r) decreases faster than 
exponentially at infinity ( [C], Lemma A.5). 

Proceeding exactly as in ||| and taking the differences 
(i) and (ii) into account, we arrive at ©• QED 

In the following, K. C R m denotes a bounded and con¬ 
nected set that is sufficiently nice so that the Poincare- 
Sobolev inequality (see |l(J, Theorem 8.12) holds on K.. 
In particular, this is the case if K, satisfies the cone prop¬ 
erty @] (e.g., if K is a ball or a cube). 

We introduce the general notation that / G L P (IC) if 

the norm \\f\\ L p(K) = [fie l/( r )N mr ] is finite. 
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LEMMA 2 (Generalized Poincare Inequality). For 
m > 2 let 1C C R m be as explained above, and let h 
be a bounded function with J^h = 1. There exists a con¬ 
stant C (depending only on 1C and h) such that for all 
sets Cl C 1C and all f £ H 1 ( 1C) (i.e., f £ L 2 (/C) and 
V/ £ L 2 (IC)) with f K fhd m r = 0, the inequality 

f a |V/(r)| 2 cTr + 7 ^ |V/(r)| 2 cTr 

>i^|/(r)| 2 cTr (12) 

holds. Here \ ■ \ is the volume of a set, and Cl c = 1C \ Cl. 

Proof. By the usual Poincare-Sobolev inequality on K 
(see m , Theorem 8.12), 

ll/lll 2 (/C) ^ C||V/|| 2 2m/( 

™+ 2 >(x) 

< 2C (|| V/|| 2 2m/ ( m+2 )(0) + 11 V/11 2 2m/(m+2 ) (QC) ) , 
if m > 2 and f K fh = 0. Applying Holder’s inequality 

||V/|| L2m/(m+2 , (n) < IIV/II^^IOI 1 /" 1 

(and the analogue with Cl replaced by Cl c ), we see that 
(|T|) holds with C = 21 K. \ 2 ' m C. QED 

The important point in Lemma 2 is that there is no 
restriction on Cl concerning regularity or connectivity. 

Proof of Theorem 1. For some R > 0 let 1C = (r £ 
R 3 , |r| < R}, and define 

</x)K = J K I^PWP Jr i I^MlVxW.fr - 

We shall use Lemma 2, with m = 3, h( r) = 
l0 GP ( r )| 2 / /r_J</’ GP | 2 , Cl = Cl^nlC and /(r) = / x (r) - 
(/x)x (see (0) and (0)). Since </> GP is bounded on JC 
above and below by some positive constants, this Lemma 
also holds (with a different constant C') with d 3 r replaced 
by |0 GP (r)| 2 d 3 r in <ED- Therefore, 


JdxJ d 3 r|^ GP (r)| 2 [/ x (r) - (/x)/cf 

< C' f dX [ f |0 GP (r)| 2 |V r /x(r)| 2 d 3 r 

/ |^ GP (r)| 2 |V r / x (r)| 2 d 3 r 
Jk 


/V -8 / 51 


. (13) 


where we used that |flx n /C| < (47 t/ 3)1V 4 / 17 . The 
first integral on the right side of © tends to zero 
as N —> oo by Lemma 1, and the second is bounded 
by ®. We conclude, since Xc |</ >GP (r)| 2 /x(r)d 3 r < 
Jr. 3 l^ GP ( r )| 2 /x(r)d 3 r, that 


iiminf <^ GP |tI^ GP > > 

AT—too N 

> f |<(> GP (r)| 2 d 3 r lim / dX [ d 3 r|T(r,X )| 2 . 
Jk n^ooJ J k 


It follows from ( |4b| ) that the right side of this inequal¬ 
ity equals [/ K |0 GP (r)| 2 d 3 r] “. Since the radius of K, was 
arbitrary, ^(^ GP |7|<(> GP ) —> 1, implying Theorem 1 (cf. 
H|, Theorem 2.20). QED 

We remark that the method presented here also works 
in the case of a 2D Bose gas. The relevant parameter to 
be kept fixed in the GP limit is g = AnN/\ ln(a 2 7V)|, all 
other considerations carry over without essential change, 
using the results in 00. A minor difference concerns 
the parameter s in Theorem 2, which can be shown to 
be always equal to 1 in 2D, i.e., the interaction energy is 
purely kinetic in the GP limit (see Jl0] ). We also point 
out that our method necessarily fails for the ID Bose gas, 
where there is presumably no BEC 0. An analogue of 
Lemma 1 cannot hold in the ID case since even a hard 
core potential with arbitrarily small range produces an 
interaction energy that is not localized on scales smaller 
than the total size of the system. There is also no GP 
limit for the one-dimensional Bose gas in the above sense. 
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